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In theories of supersymmetry breaking, it is often tlie case tliat tiiere is more tlian one metastable 
vacuum. First-order pliase transitions among such metastable vacua may generate a stochastic 
background of gravitational waves, the observation of which would provide a direct window into the 
supersymmetry-breaking sector. 



The detection of gravitational waves (GWs) will pro- 
vide unparalleled insight into the history of the early Uni- 
verse. Uninfluenced by recombination, such gravitational 
waves shed light on phenomena of the earliest epochs. In 
particular, a stochastic background of GWs is expected 
to be produced by high-energy cosmological processes, 
including first-order phase transitions. 

Such first-order phase transitions are a common pre- 
diction of new physics. Electroweak symmetry break- 
ing may entail a first-order phase transition at temper- 
atures ^ 100 GeV, while more speculative theories 
of new physics - such as grand unification or Randall- 
Sundrum models of warped extra dimensions - may pro- 
duce first-order phase transitions at temperatures rang- 
ing from T, - 1 TeV to T, - lO^^ GeV. 

Supersymmetry (SUSY), spontaneously broken at the 
electroweak scale, is among the most attractive can- 
didates for new physics beyond the Standard Model. 
The breaking of supersymmetry typically occurs in a 
sector sequestered from the Standard Model, whose 
fields inhabit a supersymmetry-breaking vacuum. This 
vacuum need not be the global minimum of the the- 
ory; indeed, metastable supersymmetry-breaking ap- 
pears to be a simple and generic feature of many SUSY- 
breaking sectors [2]. Such theories often feature one 
or more supersymmetry-breaking vacua with paramet- 
rically small rates for transition into supersymmetric 
vacua. However, it is quite possible for first-order phase 
transitions to occur rapidly among the supersymmetry- 
breaking vacua of the theory. 

In this paper, we wish to draw attention to the poten- 
tial observation of gravitational waves from phase transi- 
tions among metastable supersymmetry-breaking vacua 
in the early Universe. Such a phase transition would be 
generically expected to occur at temperatures T^, ~ \/F, 
where \/F is the primordial supersymmetry-breaking 
scale. A stochastic background of gravitational waves is 
generated when such a phase transition nears completion. 
When several bubbles of the final vacuum collide, the 
spherical symmetry of the bubbles is broken, and a frac- 
tion of the kinetic energy contained in the bubble walls is 
converted to gravitational waves. Moreover, the collision 
of two or more bubbles produces anisotropic stirring of 
the cosmic plasma, generating additional GWs. Gravita- 
tional waves from phase transitions of this genre may be 
measured at future ground- and space-based GW detec- 



tors, providing a direct window into the supersymmetry- 
breaking sector. 

Metastable Supersymmetry Breaking. ~H supersymme- 
try provides an explanation for the hierarchy problem, 
it must be broken in a sector separate from the fields 
of the Standard Model, most likely by some dynamical 
mechanism [T]. Indeed, supersymmetric gauge theories 
generically possess both supersymmetric and metastable 
supersymmetry-breaking vacua [3]. 

In a typical theory with more than one metastable 
supersymmetry-breaking vacuum, it is generally the case 
that many features of the potential separating metastable 
vacua axe set by one parameter: the primordial SUSY- 
breaking scale, \/F. In such theories, both the critical 
temperature Tc at which first-order phase transitions 
commence, and the nucleation temperature at which 
the phase transition completes, are of order ^ ^/F. Con- 
sequently, the spectrum of stochastic gravitational waves 
from a first-order phase transition among metastable 
supersymmetry-breaking vacua is characterized by VF- 

The size of the supersymmetry breaking scale \/F 
may be determined by the phenomenological require- 
ment that the MSSM soft-scale masses lie around the 
weak scale, i.e., F/M ~ 10^ - 10^ GeV, where M de- 
notes the mass scale of the mechanism mediating super- 
symmetry breaking to the MSSM. In the case of grav- 
ity mediation, M = Mp, which implies a primordial 
supersymmetry-breaking scale V^grav — 10^^ GeV. In 
the case of gauge mediation, however, the messenger scale 
M may be as low as 10^ GeV (satisfying experimen- 
tal limits on sparticle masses) and as high as 10^^ GeV 
(above which fiavor- violating contributions from gravity 
mediation become significant). Consequently, the super- 
symmetry breaking scale for gauge mediation may range 
from V^gauge — lO'^ — 10^ GeV. Within this range, low- 
scale gauge mediation (corresponding to Vf ~ 10"* — 10^ 
GeV) leads to a phenomenologically favorable spectrum 
of MSSM soft masses [4]. 

There are numerous SUSY-breaking models exhibit- 
ing more than one metastable supersymmetry-breaking 
vacuum (e.g., [3J |S]). It is often the case that the 
supersymmetry-breaking vacua of such theories are para- 
metrically close in field space, while the supersymmetric 
vacua are generated by nonperturbative effects and thus 
are sufficiently distant to make SUSY-breaking vacua 
long-lived. Moreover, for sufficiently high reheating tem- 
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peratures, thermal effects generally select vacua in which 
supersymmetry is broken [6]. 

For example, consider supersymmetric SU{Nc) QCD 
with Nf fundamental and antifundamental chiral super- 
fields Q, Q and adjoint X with polynomial superpoten- 
tial, as studied in [Xj. This theory is strongly coupled 
in the infrared, but possesses an IR-free dual provided 
Nc/2 < Nf < 2Nc/3. The adjoint scalar potential breaks 
the gauge group into a 'landscape' of vacua correspond- 
ing to the different vacuum expectation values of X. In 
the presence of a mass deformation mQQ in the super- 
potential, the low-energy theory consists of many decou- 
pled copies of the ISS model of metastable SUSY break- 
ing with gauge groups of varying size. The super- 
symmetric vacua of the theory are parametrically distant 
from the supersymmetry-breaking vacua, and tunneling 
rates into supersymmetric vacua may naturally be made 
longer than the age of the universe. The supersymmetry- 
breaking vacua, however, are separated by barriers and 
distances of order mx VF, and tunneling rates be- 
tween these vacua may be large. 

Moreover, if the Universe is reheated to a temper- 
ature TfiH ^ VF after inflation, thermal effects will 
generically select the supersymmetry-breaking vacuum 
with the largest gauge group. For the SU{Nc) SQCD 
theory with adjoints, in the infrared theory this corre- 
sponds to selecting the vacuum with the largest value 
of V^- After the metastable vacuum with the largest 
scale of supersymmetry breaking is populated, transi- 
tions to metastable vacua with smaller ^/F may rapidly 
take place [8]. As such, in these theories it is reason- 
able to expect one or more first-order phase transitions 
among different metastable SUSY-breaking vacua in the 
early universe at temperatures ~ ^/F.^ 

Gravitational Waves. ~ln general, a first-order phase 
transition will occur whenever the Universe initially finds 
itself in a false vacuum, separated from additional vacua 
of lower energy by a potential barrier. As temperature 
drops with the expansion of the universe, it becomes en- 
ergetically favorable to nucleate bubbles of a lower-energy 
final vacuum within the initial vacuum phase. Bubbles 
above a critical size expand and collide, completing the 
transition into the final vacuum. The temperature at 
which the phase transition completes may be properly 
defined as the temperature at which the probability of 
nucleating one bubble per horizon volume, per horizon 
time, approaches unity. This condition amounts to the 
requirement r(r*) = Hf, where F is the nucleation rate 
for bubbles of the final vacuum and if* is the Hubble 



parameter at temperature T*. 

As the phase transition comes to completion, collisions 
among bubbles of the final vacuum produce gravitational 
waves of characteristic frequency and fraction of the 
energy density flgw*- These gravitational waves propa- 
gate to the present era without interacting, resulting in 
an observed frequency / — f*{ai,/aQ) and observed en- 
ergy density fraction figw — ^lgw*{a^, /ao)'^{H^, / Hq)'^; a{t) 
is the cosmological scale factor. The present contribution 
of gravitational waves to the energy density as a function 
of frequency is then 
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is the energy 



where pc is the critical density and 
density in gravitational waves. 

For sufficiently strong first-order phase transitions, the 
characteristic frequency and energy density of gravita- 
tional waves rigw* produced at T* depends only weakly on 
the specific microphysics of the phase transition. Rather, 
it may be characterized by two parameters: a = e/prad, 
the ratio between the latent heat liberated in the phase 
transition and the energy density in the high energy 
phase; and /3, the rate of time variation of the nucle- 
ation rate at transition temperature T*. The strength 
of the phase transition is controlled by a, with a —^ 
and a — *■ oo corresponding to weakly and strongly first- 
order phase transitions, respectively. The duration of the 
phase transition is given by /3~^ and the size of bubbles 
by Rb ~ Wfc//?, where Vh is the velocity of the bubble wall. 

The rate per unit volume for formation of bubbles of 
the final vacuum is F ~ T'^e~^ , where S{t) is the Eu- 
clidean bounce action. Correspondingly, /3 is given in 
terms of the bounce action by /3 = H^T^ . The 

percolation condition F(r*) = Hf then implies [3] 
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Thus the quantity -ff*//3, which sets the duration of the 
phase transition and typical bubble size, is expected to 
depend only logarithmically on the scale of the phase 
transition. For the thermal effective potentials of interest 
here, a strongly first -order phase transition a ^ 1 entails 
~ Tc/Vfew, consistent with the transition tempera- 
tures implied by the percolation condition. 

Loosely speaking, the characteristic frequency of GWs 
produced by bubble collisions is given by the time scale of 
the phase transition, ~ p. Similarly, scaling arguments 
alone suggest that figw* is given in terms of /3/i?* by [TD] 



^ Whether these phase transitions give rise to topological de- 
fects depends on the particular dynamics of the supersymmetry- 
breaking sector, but the scales of interest are too low to produce 
a troublesome overdensity of defects. 



Ptot* 



{1 + ay 



T 



n 2 



(3) 



where k is the fraction of vacuum energy transformed into 
fluid kinetic energy. The spectrum of gravitational waves 
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produced by subsequent turbulent motions of the plasma 
is somewhat more complicated, and depends additionally 
on the turbulent fluid velocity Ug . In particular, the peak 
frequency of GWs produced by turbulence is set by the 
largest bubbles involved in the phase transition, and thus 
is lower than the peak frequency produced by collisions. 

For the case of detonation collisions, in which the bub- 
ble wall propagates faster than the speed of sound (as 
is the case for strongly first-order phase transitions), the 
bubble wall velocity vi,, fluid velocity Us, and energy frac- 
tion K are given simply in terms of a by I12j 
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More precise estimates of the peak frequency and en- 
ergy density stored in GWs due to bubble collisions have 
been made on the basis of numerical simulations 112j . 
More recently, an accurate analytic approximation for 
the power spectrum was obtained for gravitational waves 
produced by bubble collisions, in which the present-era 
peak amplitude and frequency of GWs from bubble col- 
lisions arc found to be [T^ 
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Here 



is the outer maximal 



fluid velocity. 

Similarly, the GW spectrum due to plasma turbulence 
has been the subject of extensive, albeit more recent, 
numerical study pPT , according to which the peak ampli- 
tude and frequency for GWs produced by plasma turbu- 
lence are given by 
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provided Us > 1/2 (as is the case for a > 0.9); for Us < 
1/2, the peak amplitude is reduced by a factor of 4^^. 

Both collisions and turbulence may contribute signif- 
icantly to the gravitational wave spectrum, with tur- 
bulence becoming dominant as the phase transition 
strengthens. The peak frequency of the resulting spec- 
trum of gravitational waves directly reflects the temper- 
ature at which the phase transition occurs, while the 
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FIG. 1: Stochastic gravitational wave signal of a first-order 
supersymmetry-breaking phase transition with = 10** 

GeV for varying strengths of the phase transition: a — 0.01 
(green), 0.1 (red), and 1.0 (blue). Dotted lines denote the 
signal from turbulence; dashed lines denote the signal from 
collisions; the solid line is the total stochastic signal for each 
value of a. 



energy density in gravitational waves depends on the 
strength of the phase transition. The above results, 
coupled with expressions for the full gravitational wave 
spectrum ilg^{f) for collisions [13^ and turbulence |14| , 
allow us to estimate GW spectra for various SUSY- 
breaking phase transitions. Characteristic gravitational 
wave spectra for first-order phase transitions with \/F ~ 
10'* GeV and various values of a are shown in Fig. [l] 

Detection. -The stochastic background of gravitational 
waves from SUSY-breaking phase transitions is likely too 
weak to be seen by ground-based bar detectors or laser in- 
terferometers such as LIGO, VIRGO, GEO and TAMA. 
On the other hand, future ground-based experiments and 
proposed space-based detectors such as LISA, BBO, and 
AGIS |15] may be sufficiently sensitive to probe phase 
transitions with a low scale of supersymmetry break- 
ing, as shown in Fig. |2] The peak sensitivity of LISA 
lies at a frequency well below the lowest frequencies ex- 
pected from SUSY-breaking phase transitions. However, 
a spaced-based atom-interferometric experiment such as 
AGIS or a correlated laser experiment such as BBO 
would be sensitive to low-scale supersymmetry breaking 
with ^/F ^ 10"* — 10^ GeV, while a future ground-based 
interferometer may probe \/F ~ 10^ — 10^ GeV. Provided 
the first-order phase transition is sufficiently strong, the 
GW signal would be visible above anticipated stochas- 
tic backgrounds from astrophysical sources and inflation. 
Unfortunately, the gravitational waves produced in theo- 
ries of gravity mediation (corresponding to peak frequen- 
cies /p ~ 1 MHz) would generally lie at frequencies too 
high to be measured by space-borne experiments. 
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FIG. 2: GW spectra alongside sensitivity curves of various 
gravitational wave detectors. From left to right, the red peaks 
correspond to GW spectra with yfF = 10*, 10^, 10*' GeV and 
Of = 1. The solid black curve denotes the sensitivity of a 
ground-based AGIS experiment (AGIS-G); the dashed curve 
denotes that of the most aggressive space-based AGIS exper- 
iment (AGIS-S). The dashed-dotted curve denotes the sensi- 
tivity of LISA, and the dotted curve the sensitivity of a corre- 
lated BBO. The AGIS sensitivity curves are taken from [15) : 
the BBO sensitivity curves from [16]. The grey band shows 
the estimated background from extragalactic white dwarf bi- 
naries |T7] , while the blue band shows the region of anticipated 
signals from inflation. 

Summary. -Stochastic gravitational waves from cosmo- 
logical processes in the early Universe provide unparal- 
leled insight into the earliest epochs. If the sector respon- 
sible for weak-scale supersymmetry breaking possesses 
more than one metastable vacuum, first-order phase tran- 
sitions among these supersymmetry-breaking vacua may 
generate a stochastic background of gravitational waves 
measurable by next-generation ground- and space-based 
interferometers. We expect a strongly first-order phase 
transition with supersymmetry-breaking scales ^/~F ^ 
10'* — 10® GeV, corresponding to low-scale gauge media- 
tion, to be accessible by future ground- and space-based 
interferometers. The measurement of such gravitational 
waves would provide a direct window into the physics of 
the supersymmetry breaking sector. 
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